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PAH emission and star formation in the host of the z~2.56 

Cloverleaf QSO 

D. Lutz 1 , E. Sturm 1 , L.J. Tacconi 1 , E. Valiante 1 , M. Schweitzer 1 H. Netzer 2 , R. Maiolino 3 , 

P. Andreani 4 , O. Shemmer 5 , S. Veilleux 6 

ABSTRACT 

We report the first detection of the 6.2/im and 7.7/im infrared 'PAH' emis- 
sion features in the spectrum of a high redshift QSO, from the Spitzer-IRS spec- 
trum of the Cloverleaf lensed QSO (H1413+117, z~2.56). The ratio of PAH 
features and rest frame far-infrared emission is the same as in lower luminosity 
star forming ultraluminous infrared galaxies and in local PG QSOs, supporting 
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a predominantly starburst nature of the Cloverleaf 's huge far-infrared luminosity 
(5.4 x 1O 12 L , corrected for lensing). The Cloverleaf 's period of dominant QSO 
activity {LboI ~ 7 x 1O 13 L ) is coincident with an intense (star formation rate 
~ 1000M o yr _1 ) and short (gas exhaustion time ~ 3 x 10 7 yr) star forming event. 

Subject headings: galaxies: active, galaxies: starburst, infrared: galaxies 



1. Introduction 



Redshifts ~2.5 witne ss both the 'quasar e poch' with peak number density of luminous 
accreting black holes (e.g. ISchmidt et al.lll995l ) and the peak in num ber of the most intens e 
star forming events as traced by the submillimeter galaxy population (jChapman et al.ll2005l ). 
suggestive of a relation of the two phenomena. Detailed evolutionary connections between 
massive starbursts and Q SOs have been discussed for many years (e.g. ISanders et al.lll988 



Norman Sz Scovilk 1988h and form an integral part of some recent models of ga laxy and 
merger evolution (e.g. lGranato et al.ll2004l ; ISpringel et al.ll2005l : iHopkins et al.ll2006l ). Phases 
of intense star formation coincident with the active phase of the quasars are a natural 
postulate of such models but have been exceedingly difficult to confirm and quantify due to 
the effects of the powerful AGN outshining tracers of star formation at most wavelengths. 

Perhaps the strongest constraint on the potential significance of star formation in QSOs 
comes from the far-infrared part of their spectral energy distribution (SEP). Indeed, this far - 



infrared emission has been interpreted as due to star formation (e.g. lRowan-Robinsonlll995l ). 
but alternative models successfully ascribe it to AGN heated dust, by postulating a dust 
distribution in which relatively cold dust at large distance from the AGN has a significant 
covering factor, for example in a warped disk configuration (ISanders et al.lll989l ). Additional 
diagnostics are needed to break this degeneracy. 



CO surveys of local QSOs (e.g. Evans et al.ll200ll : IScoville et al.ll200a Evans et alJl200q ) 
have produced a significant number of detections of molecular gas reservoirs that might power 
star formation. Depending on the adopted 'star formation efficiency' SFE=L fir/Lco the 
detected gas masses may be sufficient or not for ascribing the QSO far-infrared emission to 
star form ation. Optical studies have identified significant 'po st-starburst' stellar p opu l ations 
in QSOs JCanalizo & StocktoiJlioOll iKauffmann et al.ll2003h . On the other hand, E3 (b()05h 
suggested low star formation in QSOs, perhaps actively inhibited by the AGN, on the basis 
of observations of the [Oil] 3727A line. We have used the much less extinction sensitive mid- 
infrared PAH emission features to infer that in a sample of local (PG) QSOs, st ar formation 
is sufficient to power the observed far-infrared emission (ISchweitzer et al.ll2006l ). 
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The observational situation remains complex for high redshift QSOs. Metallicity studies 
of the broad- line region sugg est significant enrichment by star formation (e.g. lHamann fc Ferland 
19991 ; IShemmer et al.ll2004l ) but may not be representative for the host as a whole. Submm 
and mm studies of luminous radio quiet QSOs have produced significant individual detec- 
tions of d ust emission of some QSOs, as well as statistical detection of th e entire popula- 
tion (e.g. lOmont et al.ll2003l ; iPriddey et al.l 120031 ; iBarvainis fc Ivisonll2002l ). These suggest 
potential starburst luminosities up to and exceeding 10 13 L m . CO studies have detected 



large gas reservoir s in many high-z QSOs (see summaries in lSolomon fc Vanden Boutl 12005 



Greve et al.ll2005l ). Emission from high density mo l ecular gas tracers has been detected in 



some of the brightest systems (jBarvainis et al 



1997 



Solomon et al. 2003; Carilli et al. 2005 



Riechers et al.ll2006l ; iGarcfa-Burillo et al.ll2006l ; iGuelin et al.ll2007l ) and may well originate in 
dense high pressure star forming regio ns, but AGN effects on chemistry and molecular line 
excitation could also play a role (e.g. iMaloney et al.lll996l ). Finally, the [CII] 157/im rest 



wavelength fine struct ure line was detected in the z=6.42 quasar SDSS Jl 14816. 64+525150. 3 
(jMaiolino et al.l 120051 ) at a ratio to the rest frame far- infrared emission similar to the ratio 
in local ULIRGs, consistent with massive star formation. 

We have initiated a program extending the use of mid-infrared PAH emission as star 
formation tracer to high redshift QSOs. In this Letter, we use Spitzer mid-infrared spectra to 
detect and quantify star forma tion in one of the b r ightest and best stu died z~2.5 QSOs, the 
lensed Cloverleaf fH1413+117. lHazard et al1ll984J : Magain et al.lll988h . We adopt Q m = 0.3, 
Q A = 0.7 and H = 70 kms" 1 Mpc" 1 . 



Observations and Results 



We obtained low resolution (R~ 60 — 1 20) mid-infrared sp ectra of the Cloverleaf QSO 
using the Spitzer infrared spectrograph IRS (IHouck et al.ll2004l ) in staring mode on July 24, 
2006, at J2000 target position RA 14hl5m46.27s, DEC +Ild29m43.40s. The IRS aperture 
includes all lensed images. 30 cycles of 120sec integration time per nod position were taken 
in the LL1 (19.5 to 38.0 /im) and 15 cycles in the LL2 (14.0 to 21.3 jttm) module, leading to 
effective on-source integration times of 2 and 1 hours, respectively. We use the pipeline 14.4.0 
processed basic cali brated data, own deglitching and coaddition procedures, and SMART 
( jHigdon et al.ll2004j ) for extraction. When combining the two orders into the final spectrum, 
we scaled the LL2 spectrum by a factor 1.02 for best match in the overlapping region. 

Fig. [1] shows the IRS spectrum embedded into the infrared to radio SED of the Clover- 
leaf, and Fig. [2] the IRS spectrum proper, together with the location of key features in the 
corresponding rest wavelength range. The rest frame mid-infrared emission is dominated by 
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a strong continuum, approximately flat in uF v , due to dust heated by the powerful active 
nucleus to temperatures well above those reached in star forming regions. Superposed on 
this continuum are emission features, which we identify with the 6.2/zm and 7.7/xm aromatic 
'PAH' emission features normally detected in star-forming galaxies over a very wide range 
of properties. As expected for a Type 1 AGN, there are no indications for the ice (6/im) 
or silicate (9. 6/im) absorptions seen in heavily obscured galaxies. None of the well-known 
emission lines in this wavelength range is bright enough to be significantly detected in this 
low resolution spectrum, although we cannot exclude a contribution of [Ne Vh 1 7.64/im to 



the 7 . 7nm feature. Adopt ing standard mid-infrared low resolution diagnostics (IGenzel et al. 



19981 ; lLaurent et al.ll2000l ). the weak PAH features on top of a strong continuum agree with 
the notion that the Cloverleaf is energetically dominated by its AGN. The detection of PAH 
features with several mJy peak flux density in a z~2.6 galaxy, however, implies intense star 
formation, which we discuss in conjunction with other properties of the Cloverleaf. 

By fitting a Lorentzian superposed on a local polynomial continuum, we measure a 
flux of 1.52xlO~ 21 Wcm~ 2 for the 6.2/im feature, with a S/N of 6. The 7.7//m feature is 
more difficult to quantify. ISchweitzer et al.l (120061 ) have discussed PAHs as star formation 
indicators in loc al (PG) QSOs, PAHs are also detected in the average QSO spectrum of 
Hao et al.l (120071 ). The AGN continuum of those QSO s shows superpos ed silicate emission 
features at > 9/zm (see also ISiebenmorgen et al.l 120051 ; lHao et al.l 120051 ). If PAH emission 
is additionally present, the PAH features partly 'fill in' the minimum in the A GN emission 
before the onset of the silicate feature (see Fig. 2 of ISchweitzer et al.l (120061 )). causing a 
seemingly flat overall spectrum. In reality, there is simultaneous presence of AGN continuum, 
the 6.2-8.6/xm PAH complex, silicate emission, and more PAH emission at longer wavelengths. 
From inspection of Fig. [2], a similar co-presence of PAH and silicate emission is observed for 
the Cloverleaf. We note that the presence of silicate emission in the luminous Cloverleaf 
Type 1 QSO, o ther high z Type 1 QSOs (Maiolino e t al., in prep) as well as in luminous 
Type 2 QSOs Jsturm et all bood : iTeplitz et all IgOOgh has implications for the location of 



Hao et al 



20051 ) in unified AGN schemes. We leave 



this cool silicate component (~200K, 
further discussion of the properties of silicates to a future paper. For the 7.7/im feature of 
the Cloverleaf, we adopt a flux of 6.1x lO~ 21 Wcm~ 2 . Th is flux was determined by scaling a 
PAH template (ISO spectrum of M82. ISturm et al.ll2000l ) to the measured Cloverleaf 6.2/im 
feature flux, and then fitting three lorentzians to represent the 6.2, 7.7, and 8.6 features plus 
a local polynomial continuum to this tem plate. Similar Loren tzian fits were also used for 
local comparison objects discussed below. iBrandl et al.l (120061 ) have quantified the scatter 
of the 6.2 to 7.7/iin flux ratio in starbursts, with 0.07 in the log this scatter indicates the 
modest uncertainty induced by tying the longer wavelength features to the 6.2/xm one. The 
result of subtracting the scaled M82 template from the Cloverleaf spectrum is indicated 
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in Fig. [21 and shows a combination of continuum and silicate emission very similar to local 
QSOs. Directly measuring the 7.7/zm flux by fitting a single Lorentzian plus local continuum 
to the Cloverleaf spectrum gives a ~40% lower feature flux, which would be a systematic 
underestimate because of the complexity of the underlying continuum/silicates discussed 
above. 



3. Intense star formation in the host of the Cloverleaf QSO 



The Cloverleaf SED (Fig. [T]) shows strong rest frame far-inf rared emission in addition 
to the AGN heated dust emitting in the rest frame mid-infrared. IWeiB et al.l (120031 ) decom- 
posed the SED into two modified blackbodies of temperature 50 and 115 K, the rest frame 
far-infrared (40-120/xm) luminosity of 5.4 x 1O 12 L is dominated by the colder component 
and could largely originate in star formation. Comparison of PAH and far-infrared emission 
can shed new light on this question. The bolometric (rather than rest-frame far-infrared) 
luminosity of the Cloverleaf will still be dominated by the AGN. We estimate Lb i extrap- 
olating from the observed rest frame 6/im continuum which f or a mid-infrared sp ectrum 
with weak PAH but strong continu um will be AGN do minated ( Laurent et al.ll2000l ). Using 



LboI ~ 10 x uL u (6fim) based on an lElvis et al.l (119941 ) radio-quiet QSO SED, the AGN lu- 
minosity is ~ 7 x 1O 13 L . A si milar estimate ~ 5 x lO 13 !^ is obtained from the rest frame 



optical (observed near-infrared; iBarvainis et al.l 119951 ) continuum, tracing the AGN ionizing 
continuum, and the same global SED. 



Schweitzer et al.l (120061 ) have measured PAH emission in local QSOs and compared the 



PAH to far-infrared emission ratio to that for starbursting ULIRGs, i.e. those among a larger 
ULIRG sample not showing evidence for dominant AGN and not having absorption domi- 
nated mid-infrared spectra. Fig. [3]places the Cloverleaf on their relation between 7.7/im PAH 
luminosity and far-infrared luminosity. L(PAH)/L(FIR) is 0.014 for the Cloverleaf, very close 
to the mean value for the 12 starburst-dominated ULIRGs of < L(PAH)/L(FIR) >= 0.0130. 
The scatter of this relation is 0.2 in the log for these 12 comparison ULIRGs, indicating the 
minimum uncertainty of extrapolating from the PAH to far-infrared emission. The Clover- 
leaf thus extends the relation between PAH and far-infrared luminosity for the local QSOs 
and ULIRGs to ~5 times larger luminosities. Its PAH emission is consistent with an ex- 
tremely luminous starburst of ULIRG-like physical conditions powering essentially all of the 
rest frame far-infrared emission. 



Teplitz et al.l (120061 ) present the IRS spectrum of the lensed FIR-bright Type 2 AGN 



IRAS F10214+4724 at similar redshift. They report a marginal feature at 6.2/zm rest wave- 
length which they do not interpret as PAH given the lack of a 7.7/mi maximum. Given 
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the strength of silicate emission in this target, PAH emission may be present in the blue 
wing of the silicate feature without producing a maximum, and such a component may be 
suggested by comparing their Fig. 1 with the later onset of silicate emission in the spectra 
of local QSOs. The tentative 6.2/im peak in IRAS F10214+4724 has similar peak height as 
the Cloverleaf PAH feature, in line with our interpretation and the similar rest frame FIR 
fluxes of the two objects. 

With ~5-10 % of its total luminosity originating in the rest frame far-infrared and by star 
formation, the Cloverleaf is within the range of local QSOs, and not a pronounced infrared 
excess object. Specifically, its ratio of FIR to total luminosity and the ratio of res t frame far- 
i nfrar ed (60/xm) to mid-infrared (6/xm) continuum are ab out twice those o f the lElvis et al. 
(119941 ) radio-quiet QSO SED. Adopting the conclusion of ISchweitzer et al. that star 



formation already dominates the FIR emission o f local PG QSOs an d considering the modest 
FIR 'excess' of the Cloverleaf compared to the lElvis et al.l (119941 ) SED then suggests only 
a small AGN contribution to its FIR luminosity. Other z~2 QSOs may have lower ratios 
of FIR and total luminosity, and conversely larger AGN contributions to their more modest 
FIR emission, though. After correctin g for lensing, the Clov erleaf submm flux is a factor ~2 
above the typical bright z~2 QSOs of iPriddey et al.l (120031 ) whose rest frame B magnitudes 
in addition are typically brighter than the delensed Cloverleaf. 

Sub millimeter galaxie s host starbursts of s i milar luminosity as the Cloverleaf, at similar 
redshift. iLutz et al.l (120051 ) and lValiante et al.l (120071 ) have obtained IRS mid-infrared spec- 
tra of 13 SMGs with median redshift 2.8, finding mostly starburst dominated systems. A 
comparison can be made between PAH peak flux density and flux density at rest wavelength 
222yum which is obtained with minimal extrapolation from observed SCUBA 850/im fluxes. 
Combi ning the 7. 7/im feature pea k of 5.1mJy (Fig. [2]) with a z/ 3 5 extrapolation of the SCUBA 
flux of iBarvainis fc Ivisonl (120021 ) places the Cloverleaf at Log(Si MiiLz/S222jxm3 ~ — 1-2, near 



the center of the distribution of this quantity for the SMGs of IValiante et al.l (120071 . their 
Fig. 4). Like the SMGs, the Cloverleaf appears to host a scaled up ULIRG-like starburst, 
but with superposition of a much more powerful AGN, also in comparison to the gas mass. 

Tracers of high density gas, in particular HC N but also HCO + have been detected in a 



few high redshift QS Os including the Cloverleaf (iBarvainis et al.lll997l ; ISolomon et al.l 12003 



Riechers et al.l 120061 ). Their ratio to far- infrared emission is similar to the one for Galactic 
dense star forming regions, and has been used to argue for dense, high pressure star forming 
regions dominating t he far-infrared luminosity of these QSOs as well as of local ULIRGs (e.g. 
Solomon et al.l 120031 ) . Intense HCN emis sion is observed also from X-ray dominated regions 



close to AGN (e.g. iTacconi et al.l Il994l ). and there is ongoing debate as to the possible 



contributions of chemistry and excitation in X-ray dominated regions, and other effects like 
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radiatiye pumping, to the emission of dense molecular gas tracers in ULIRGs and QSOs 
jKohnolliooa llmanishi et alibood : iGracfa-Carpio et alJlioOq l. Unlike HC N, PAH emission 
is severely reduced in X-ray dominated regions close to AGN (jVoitl Il992l ) and provides an 
independent check of the effects of the AGN on the molecular gas versus the role of the host 
and its star formation. In a scenario where XDRs dominate the strong HCN emission and 
the hosts PAH emission, reproducing the consistent ratios of these quantities to rest-frame 
far-infrared over a wide range of far-infrared luminosities would thus require a considerable 
amount of fmetuning. In contrast, these consistent ratios are a natural implication if all 
these components are dominated by ULIRG-like dense star formation. 

Our detection of PAH emission is strong support to a scenario i n whi ch the Cloverleaf 
QSO coexists with intense star formation. Applying the iKennicuttl (119981 ) conversion from 
infrared luminosity to star formation rate to Lfir = 5.4 x 1O 12 L suggests a star formation 
rate close to 1000 M yr _1 , which can be mainta ined for a gas exh austion timescale of only 
3 x 10 7 yr, for the molecular gas mass inferred by lWeifi et al.l (120031 ). At this time resolution, 
the period of QSO activity coincides with what likely is the most significant star forming 
event in the history of the Cloverleaf host. 



This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology, under a con- 
tract with NASA. Support for this work was provided by NASA under contracts 1287653 
and 1287740 (S.V.,O.S.). We thank the referee for helpful comments. 
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Table 1. Clover leaf properties 



Quantity 


Value 


Reference 


Redshift z 


2.55784 


Weifi et al. (20031 


Amplification /i^ 


11 


Venturini & Solomon (20031 


F(PAH 6.2/xm) 


1.5 x lCr 21 Wcm- 2 


this work 


F(PAH 7.7/im) 


6.1 x 10- 21 Wcm- 2 


this work 


L(PAH 7.7/im) a 


7.6 x 1O 1O L 


this work 


L(40-120/mi) a 


5.4 x 1O 12 L 


WeiB et al. (2003) 


M(H 2 ) a 


3.0 x 1O 1O M 


WeiB et al. (2003) 


WQSO) a 


~ 7 x 1O 13 L 


this work, 10 x uL v (6fim) 



a Corrected for lensing amplification 11 and to our adopted cosmology Q m = 0.3, Q\ = 0.7 
and H = 70 kms^ 1 Mpc" 1 (D L =20.96 Gpc). 



- 12 - 



1.0000 



0.1000 - 



i 0.0100 

"D 

X 

3 



0.0010 - 



0.0001 



10 L 
— r 



Rest wovelength [yum] 



10 

— r 



10 z 
— r 




_i I I 



_i I I 



_i I I 



10' 
— r 



10" 
— r 



10' 



Cloverleaf 
(H1413+1 17) 



1 



_i I I 



10 u 10 10^ 10° 10 

Observed wavelength [/xm] 



10' 



10" 



Fig. 1. — Infrared to radio spectral energy distribution for the Cloverleaf QSO. The 



(Barvainis et al. 
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flux appears too high while the other mid-infrared data are consistent within plausible cali- 
bration uncertainties. 
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Fig. 2. — IRS spectrum of the Cloverleaf QSO. The PAH emission features as well as the 
expected locations of strong spectral lines in this wavelength range are marked. The dotted 
line shows the spectrum after the subtraction of a PAH template (spectrum of M82, see also 
bottom of figure), redshifted and scaled to the measured strength of the Cloverleaf 6.2/iin 
PAH feature. Note that the noise in IRS low resolution spectra increases strongly from 
~33/iin towards the long wavelength end. 
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Fig. 3. — Relation of 7.7/im PAH luminosity and rest fra me FIR luminosity for t he Cloverleaf 
and for local PG QSOs and starbursting ULIRGs from ISchweitzer et al.l (12006I ). 



